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Abstract 

Consider a spherically symmetric spacelike slice through a spherically sym- 
metric spacetime. One can derive a universal bound for the optical scalars 
on any such slice. The only requirement is that the matter sources satisfy 
the dominant energy condition and that the slice be asymptotically flat and 
regular at the origin. This bound can be used to derive new conditions for 
the formation of apparent horizons. The bounds hold even when the matter 
has a distribution on a shell or blows up at the origin so as to give a conical 
singularity. 
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Relativists, especially those who are numerically inclined jlj , have long known that reg- 
ular spacelike slices often wrap around singularities rather than approaching them. In this 
letter we derive a new and remarkable relation which gives a bound on the optical scalars 
and which shows how slices which are asymptotically flat may be prevented from coming 
close to singularities. 

Consider a spacelike slice through spacetime. The geometry of this slice cannot be 
chosen at will; it must satisfy the constraint equations. In the spherically symmetric case 
these constraints can be written as equations for the optical scalars. These equations, 
combined with the requirement of regularity at the origin and at infinity, force the optical 
scalars to remain bounded over the entire slice. The optical scalars are four- dimensional 
objects which we expect to become unboundedly large as one approaches a singularity. 
Thus regular spacelike slices are excluded from regions near singularities. This bound on 
the optical scalars also has a more immediate use. Over the years, we have been interested 
in developing criteria to determine when and if apparent horizons form [0|§. In spherically 
symmetric systems the existence of an apparent horizon implies the existence of a black hole 
^H]. These bounds on the optical scalars allow us sharpen significantly our condition for 
the formation of apparent horizons. 

We define a spherically symmetric spacetime as one having the metric 

ds 2 = -a 2 (r, t)dt 2 + a(r, t)dr 2 + b(r, t)r 2 dtt 2 (1) 

where < < 2tt and < 9 < n are the standard angle variables such that dVt 2 = 
d6 2 + sin 2 6d<p 2 . 

The initial data for the Einstein equations are prescribed by giving the spatial geometry 
at t — 0, i.e., by specifying the functions a(r, 0) and 6(r, 0), and by giving the extrinsic 
curvature (again at t — 0) 

K = = «-=M^l (2) 

2aa v an ^Jaba 

where the areal (Schwarzschild) radius R is defined as 



R = rVb. (3) 

It is useful to define the mean curvature of a centered sphere in the initial hypersuface by 

2d r R 

P = ^R- () 

In a general spacetime the behaviour of a pencil of lightrays is described by specifying 
a number of functions which describe the expansion and shear of the rays. In a spherically 
symmetric spacetime, however, we need specify only two. These objects can be expressed in 
terms of the initial data on any spacelike slice so they are simultaneously three-dimensional 
and four- dimensional scalars. These optical scalars are the divergence of future directed 
light rays 

B = W R = p-K; + trK, (5) 
K aat out 

and the divergence of past directed light rays 

e ' = -W^r + R = p + K-trK, (6) 
K aatin 

el £i £b i-l fi si 

where —77. = — — ■%= and —7- = — + ■% are the full derivatives in the direction orthogonal 

aat in a v a aat out a v a 

to the centered sphere of ingoing and outgoing photons respectively. In flat space-time both 
quantities are positive and equal to 2/R, where R is the radius of a sphere; hence each of 
the products R9 and R9' equals 2. 

The initial data must satisfy the constraints. These constraints, expressed in terms of 6 
and 9', can be written as 

d^R) = -8nR(p - i?=) - -^ E (9 2 R 2 - 4 - A9trKR 2 
y/ a AR ^ 

+9R(9R-9'R)), (7) 

d l {9'R) = -8tt R(p + 4=) - -^{0' 2 R 2 - 4 + 49'trKR 2 

y/ a 4:R ^ 

+9'R{9'R-9R)), (8) 
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where I is the proper distance from the center, i.e., dl = ^Jadr. p and j r are the energy 
density and the current density of the sources that generate the gravitational field. Note 
that jr/y/a equals j.n where n is the unit normal in the radial direction. We will assume 
that the sources satisfy the dominant energy condition, p > If the origin is regular, local 
flatness forces both optical scalars to satisfy the conditions lim^o OR = lim^o O'R = 2. 
Asymptotic flatness also gives lim^oo 9R = lim^oo 9'R = 2. 

The primary result of this calculation is a proof that if OR, Q'R are bounded at the origin 
and at infinity they are bounded on the entire hypersurface. Define B = 4 sup 0<r<oo (|i?tr_ft'|). 
We prove 

Lemma 1. Given spherical initial data that are regular at the origin and at infinity 
with sources that satisfy the dominant energy condition, both optical scalars are bounded 
on the entire hypersurface 

- 2 - B < OR, 9'R<2 + B. (9) 

Proof: Let us assume that OR > 2 + B and OR > 9'R. Consider the non-source part of 
eqn.(0), i.e., (0 2 R 2 - 4 - AOtrKR 2 + 8R(8R - O'R). Since 0R>2 + B, the first three terms 
are nonnegative while OR > O'R means that the last term is nonnegative. Therefore eqn.(^) 
implies that di(0R) < 0. Also, if O'R > 2 + B and O'R > OR, a similar analysis of eqn.(||) 
gives di{0'R) < 0. Therefore, if max(^i?, 9'R) >2 + B then di[m&x(9R, 9'R)] < 0. Since the 
maximum starts at 2, and the derivative at 2 + B is negative, it cannot rise above 2 + B. 
Hence 2 + B is an upper bound. 

The argument for the lower bound works in exactly the same way. Let us assume that 
9R < -2 - B and that 9R < 9'R. Again eqn.(0) means that d t {0R) < 0. Hence if 
min(9R,9'R) < —2 — B then di[mm(9R,9'R)] < 0. This means that one or the other of 
(9R, 9'R) is driven more and more negative. However, asymptotic flatness demands that 
both rise up to +2 at infinity. Contradiction! 

Let us stress that while B is a three-scalar which depends on the particular spacelike slice, 
9, 9' and R are four-dimensional scalars, properties of the spacetime, which are independent 
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of the choice of foliation. Thus eqn.([|) places restrictions on the kind of regular spacelike 
slice that may enter particular regions of spacetime. 

There are only two allowed topologies for globally regular, asymptotically flat, spherically 
symmetric, spacelike three-manifolds. They can either have R 3 topology with a regular 
center and one asymptotic end or R x S 2 topology with two asymptotic ends, as in the 
Schwarzschild geometry. Lemma 1 holds in both cases. 

Lemma 1 has a number of interesting consequences. Let us assume, for a moment, that 
the trace of the extrinsic curvature vanishes, i.e., that the initial data define a maximal 
slice. This means that B = and Lemma 1 implies that \6R\, \6'R\ < 2. A surface on which 
6 < is called a trapped surface; such surfaces play a key role in the singularity theorems 
of general relativity. Eqn.(^) can be used to derive 

di{9R 2 ) = -8nR 2 (p - A) + i + Ur(26'R - OR). (10) 

V CL 4 

Let L(S) be the geodesic (proper) radius of a sphere S; R(S) its areal radius; M(S) = 
Iv(s) P&V the total mass inside S; and P(S) = fv(S) ^a^V be the total radial momentum. 
Integrating ( [TDD noting that AiiR 2 dl = AiryfaR 2 dr is the proper volume, we get 

{6R 2 ){S) = -2{M - P){S) + L(S) 
1 rHS) 

+- / 9R(26'R-9R)dl. (11) 
4 Jo 

We can see that J J L dWR{26'R -9R)<\ J L dl(6'R) 2 < L, where the first inequality comes 
from the trivial estimate 2ab — a 2 < b 2 and the second from Lemma 1. Therefore 

(0R 2 ){S) < -2{M-P){S) +2L(S) (12) 

for any surface S. In particular, if M — P > L at any given sphere S then 0(S) must be 
negative. Thus we have proven: 

Theorem 1. Under conditions of Lemma 1, assuming trK = 0, if the difference between 
the total rest mass M(S) and the radial momentum P(S) exceeds the proper radius L(S) 
of a sphere S, M(S) — P(S) > L(S), then S is trapped. 
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This theorem improves our earlier result 0, in which we got a similar result but with 
L replaced by |L and the weaker conclusion that there exists a trapped surface inside S. 
The difference is due to the fact that we now impose the somewhat stronger condition that 
P — \ j\ > 0, whereas in || we used p + ^(K^) 2 > 0. Since the new conditions in Theorem 
1 eliminate tachyons this is a real difference. The constant | also appears in our criteria for 
the formation of cosmological black holes ||; we believe that these can also be improved to 
1. 

The meaning of Theorem 1 is transparent. Radially ingoing matter j r < helps form 
apparent horizons. The presence of outgoing matter, i.e ., when P(S) becomes positive, 
has to be compensated for by a greater matter density. In the extremal case of radially 
outgoing photons, when M(S) = P(S), apparent horizons cannot form. This follows from 
our Theorem 2 below. 

Theorem 1 is sharp in the sense that there exists an initial value configuration when 
the inequality saturates. This is a 3-geometry created by a shell of moving matter; the 
explicit calculation will be done elsewhere. The case in which P = was discussed in |2j 
and the corresponding criterion (with the same constant 1, as above) was shown to be the 
best possible. 

It is interesting that we obtain an exact criterion with the constant 1; this suggests that 
Theorem 1 is part of a more complex true statement that can be formulated for general 
nonspherical spacetimes. It suggests also that M(S) is a sensible measure of the energy 
of a gravitational system that might appear as a part of a quasilocal energy measure in 
nonspherical systems. 

We also obtain a necessary condition for the formation of apparent horizons. In || 
we found a criterion based on asymptotic data outside a collapsing system, [fj] states that 
M(S) > ~ must be satisfied if S is trapped in the case of moment of time symmetry data. 
The same holds true if the matter is moving under some stringent conditions on the sign of 
the momentum density Here we will derive a different (and not particularly interesting, 
although exact) estimate. The most important assumption we make is that 9' is everywhere 
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positive on the initial hypersurface. Just as 9 < guarantees a singularity to the future, 
9' < guarantees a singularity to the past. Therefore, data which arises from a regular past 
must have positive 9'. 

Theorem 2. Assume a regular maximal slice on which the sources satisfy the dominant 
energy condition. Let S be the innermost trapped surface and let (R9') > e > inside S. 
Then 

M{S) - P{S) > "-L. 
Proof. As before, we consider ([0]), which reads 

9R 2 = -2(M -P)+L+- [ L dWR(29'R - 9R). (13) 

4 jo 

Inside S, R9 is positive. We seek a lower bound on the last term on the right hand side of 



(|13|). Let t = R9,u = R9'; from Lemma 1 we know | t |, | u |< 2, so our task consists in 
estimating 2tu - t 2 for < t < 2, e < u < 2. We know that 2tu - t 2 > F(t) = 2te - t 2 . The 
only extremum of F(t) is a maximum at t — e. The minimum must occur at the endpoints 
and it is easy to show that 2tu — t 2 > F(t) > 4e — 4. Inserting this into (|13D yields 

9(S)R 2 >-2(M- P) (S) + L(S) + - / dlUe - 4) 

4 Jo 

= -2(M - P)(S) + eL, (14) 

that is, since 9(S) = 0, 

M(S)-P(S)>f. 

Hence Theorem 2 is proven. 

The inequality of Theorem 2 becomes an equality in the case of a spherical shell. The 
geometry inside the shell is flat and 9'R = 2. The necessary condition that the shell be 
trapped is that M — P > L. In we proved this in the special case when P = 0. 

It is clear that the analysis performed here can include cases where the sources are 
distributions rather than classical functions; in particular, we have no difficulty with shells 
of matter. All we get on crossing the shell is a downward step in 9 and 9' . More interestingly, 
we can extend the analysis to include weak singularities at the origin. 



Let us begin by considering a conical singularity ||. Consider a metric of the form 



dS 2 = dr 2 + a 2 r 2 dQ 2 . 



(15) 



The scalar curvature of this metric is = 2(1 — a 2 )/r 2 . A moment of time symmetry data 
set is one for which j l and = 0. For such data sets the constraints reduce to ( 3 >R = Wnp. 



For the above metric we get p = (1 — a 2 )/87rr 2 . The dominant energy condition reduces to 
the positivity of p, which implies a 2 < 1. For this metric we can also compute the mean 
curvature p, which in this case equals both 9 and 6', to get p = 2/r = 2a/R. Hence we get 
\pR\ < 2. However, the argument of Lemma 1 only requires that 8R, 8'R be bounded at the 
origin. Therefore we have shown that Lemma 1 holds for moment of time symmetry data 
with a conical singularity at the origin. The conical singularity in question is determined 
by the deficit of the solid angle 4tt(1 — a 2 ). We will show that a similar result holds true for 
general nonmaximal data. 

Let us consider initial data such that trK is finite while R9 — > X and R9 1 — > Y as 
R — > 0. Let us also assume that di(R8) and di(RO') are finite at R — 0. There are terms 
on the right hand side of eqns.(7) and (8) which seem to diverge like 1/R. The source term 
will have the same sort of 1/R divergence if 8irR 2 p — > a and 8iiR 2 j r / yfa — > ft, just as in the 
case of the conical singularity. The coefficient of this 1/R term must vanish. This gives us 
a pair of equations, one from (7) and one from (8) 



a - (5 + X -X 2 




1 



0; 



(16) 



a + (3 + -Y 2 - -XY 
2 4 



1 



0. 



(17) 



By adding these equations we get 



4a = 4 - X 2 - Y 2 + XY; 



(18) 



and by subtracting 



A(3 = X 2 - Y 2 . 



(19) 



S 



Note that eqn. fllED implies that a < 1. The weak energy condition gives a > \(3\. Let 
us assume that (3 > 0. Eqn. ([T9|) now gives us X 2 > Y 2 and Y = ±\/X 2 — 4(3. Substituted 
this into eqn.flTBD to give 



[3X 2 - 4(1 - a + (3)} [X 2 - 4(1 - a + (3)} + AX 2 (3 2 = 0. (20) 

The roots of this equation, if it has any, must lie in the range 4(1 — a + j3)/3 < X 2 < 
4(1 — a + (3). Therefore we have shown that 2 > \X\ > \Y\. If we assume (3 < 0, we just 
reverse the roles of X and Y. Hence we obtain 

Lemma 2. Given p > \j\ and if all of trK, OR, 6'R, d t 6R, diO'R, (8vr f* P R 2 dR)/R are 
finite in the limit R = then 

2 > lim \9R\, Hm \9'R\, x> ^^pR 2 dR { ) 

From eqns.(4), (5) and (6) it is clear that 

„ 2d r R n OR + 6'R 

2d l R=^ r =pR= . 22 

\/a 2 



This means that the spatial part of the metric (1) can be written, at least in a small 
neighbourhood of R = 0, as 

wfm dRl+R2ditK (23) 

The estimate derived in lemma 2 implies that, under the stated conditions, there can be 
at most a conical singularity at the origin, with solid angle deficit 47r(l — ^ X ~^p )■ Conical 
singularities have previously been investigated in 2+1 gravity 0. In the 2+1 case the 
conical singularity can also be described by an angle deficit expressed in terms of the mean 
curvature: 27r(l —pR). However, in the 2+1 case the geometry is locally flat but globally 
nontrivial and the deficit angle is related to a total mass |J. In our case, the deficit angle 
is a local phenomenon caused by a mildly singular mass distribution at the origin, where p 
diverges like r~ 2 . 

Lemma 2 gives the desired bound \9R\, \9'R\ < 2 at the origin so we get a generalized 
version of Lemma 1: 
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Lemma 1'. Assume an asymptotically flat nonmaximal slice, satisfying the dominant 
energy condition, such that 4sup 0<JKoo \RtrK\ = B is finite. Let the conditions of Lemma 
2 be satisfied at the origin. Then 

2 + B>\6R\,\6'R\. (24) 

Theorems 1 and 2 hold under similar conditions. 

As we have mentioned earlier, OR and 6'R are defined for any point in a spherically sym- 
metric spacetime geometry, independent of any foliation or choice of time. One consequence 
of Lemma 1 is that if a point exists in a spherical spacetime for which either \0R\ or \0'R\ 
is larger than 2 then we know that a regular, maximal, asymptotically flat slice cannot pass 
through this point. 

Consider regular, asymptotically flat, spherically symmetric initial data which contain 
an apparent horizon. Let us now evolve the spacetime and look at the maximal Cauchy 
development of this data. We are guaranteed that a singularity will occur for a sufficiently 
large value of local proper time. It seems to us that there are only three realistic outcomes: 

i) The singularity will be of the kind where R9 — > — oo, as in the Schwarzschild singularity. 
Maximal slices (and any other slicing with a regular trace of the extrinsic curvature) do not 
cover the full Cauchy evolution. We get a "collapse of the lapse" . The foliation can continue 
for infinite time as seen by asymptotic observers but "freezes" near the interior. 

ii) We get some sort of "bag of gold" forming, where R goes to zero at some finite value 
of the proper radius, L, and part of the spacetime pinches off from the rest. 

iii) A singularity appears with diverging mass density. This may be a shell-crossing 
singularity, a central conical singularity as we discussed above or some sort of 'strong' central 
singularity. We expect that the appearance (or otherwise) of these singularities would be 
determined by the matter equation of state. 
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